Odd-even binding effect from random two-body interactions 
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Systematic odd-even binding energy differences in finite metallic particles are usually attributed 
to mean-field orbital energy effects or to a coherent pairing interaction. We show analytically and 
numerically that a purely random two-body Hamiltonian can also give rise to an odd-even staggering. 
We explore the characteristics of this chaotic mechanism and discuss distinguishing features with 
respect to the other causes of staggering. In particular, randomness-induced staggering is found to 
be a smooth function of particle number, and the mechanism is seen to be largely insensitive to the 
presence of a magnetic field. 
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I. INTRODUCTION 

Interacting finite fcrmionic systems such as atomic nu- 
clei, metallic clusters and small metallic grains j^] 
display an odd-even staggering in ground-state energies, 
i.e., the binding energy of an even- number system is 
larger than the arithmetic mean of its odd- number neigh- 
bors. There are two well-known mechanisms that can 
give rise to this staggering, namely the Kramers degen- 
eracy in the mean- field Hamiltonian and the BCS mech- 
anism arising from an attractive effective interaction. In 
nuclei, the BCS pairing mechanism resulting from a resid- 
ual nucleon-nucleon interaction is dominant but the 
mean-field or orbital energy effect may also be signif- 
icant in the lighter nuclei [Q. Surprisingly, many ba- 
sic phenemona normally associated with pairing can also 
arise from random interactions. The behavior of random- 
interaction ensembles has mostly been studied in a nu- 
clear physics context |p|-pl| but there has also been some 
study of spectra in the context of small metallic grains 
& 

In the case of metallic clusters of fewer than a hun- 
dred atoms, the orbital energy effect is rather strong and 
staggering is seen for species that do not exhibit super- 
conductivity. This effect can be easily understood using 
a jellium model or density functional theory |13]. On the 
other hand, the staggering effect seen in Ref. [ |] may have 
some contribution from the BCS pairing mechanism. A 
number of theoretical studies have been made [M us- 
ing techniques applicable to large finite systems [fl5|| . 
Taking a uniform mean-field spectrum and an attractive 
pairing interaction with constant coupling, one observes 
a smooth crossover from BCS superconductivity in the 
bulk to the few-electron regime. For small systems, the 
gap is of the size of the mean level spacing and thus ceases 
to be an indicator for pairing. Nevertheless, strong pair- 
ing correlations and odd-even staggering persist as the 
system size decreases. 

In a grain with irregular boundaries, one expects that 



the electron orbitals will have a chaotic character and 
therefore the interaction will have a random as well as 
a regular part. In this paper we will introduce such an 
interaction and study its typical effects on the binding 
systematics. Our Hamiltonian thus includes attractive 
and repulsive pairing interactions as well as more general 
two-body interactions. The assumption of randomness 
is motivated as follows: For nuclei it is well known that 
the residual interaction leads to fluctuation properties in 
wave functions and energy levels that are similar to those 
of random matrices taken from the Gaussian orthogonal 
ensemble [H. In the case of small metallic grains or 
quantum dots, one may assume that their irregular shape 
leads to chaoticity in the single-particle wave functions 
|]l7| . This in turn causes randomness in those two-body 
matrix elements that link four different orbitals with each 
other. Matrix elements between pairs of orbitals that are 
related by time-reversal symmetry need not necessarily 
be random, and these determine the "coherent" terms of 
the interaction. 

A realistic Hamiltonian for quantum dots or small 
metallic grains would thus conserve total spin and in- 
clude spin-independent one-body terms, random two- 
body interactions, and coherent interactions that are 
non-random but have attractive and repulsive compo- 
nents. The most general Hamiltonian to study generic 
properties when all these features are included may be 
written as 
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+vi, ijk i(<Ji\a\(T2) ■ (<J3\a\a 4 )}cl ai c] a2 c ka3 ci a4 . (1) 

Here the coherent parts of the interaction are repre- 
sented by the terms with coefficients u, w s , and g. The 
fluctuating parts of the interaction are represented by 
the terms containing u' , w' s , g', and v s . These fluctuating 
parts are typically taken from ensembles with a Gaussian 
distribution; they are thus characterized by the width of 
the Gaussian. The single-particle term sets the energy 
scale and may often be taken to give a uniform spacing of 
levels without loss of generality. This full Hamiltonian is 
difficult to study due to its many parameters. There have 
been many studies in the limit in which fluctuation ef- 
fects are only included in the single-particle Hamiltonian 
£i ]l8|-pH. We consider a very different limit, neglecting 
the coherent terms in the interaction and assuming the 
v s term to dominate the fluctuating parts. Properties of 
such random two-body interaction ensembles have been 
studied extensively in nuclear physics |22]-|2"5|] . 

When the Hamiltonian of the nuclear shell model was 
modeled in this way, it was found that the spectral prop- 
erties were quite regular for the ground states. As ex- 
amples we mention J p = + ground-state dominance in 
shell model calculations with random interactions || ||, 
band structure in interacting boson models with random 
couplings structure in ground-state wave functions 
of two-body random ensembles ]2(| , and an odd-even 
binding effect in filling a large shell In the context 
of quantum dots, the random two-body interactions were 
found strongly to favor singlet ground-state spins Jl2|,p7t ■ 
Recently, this structure has been investigated using the 
group symmetry of the random Hamiltonians p8[ | . These 
findings suggest that the structure of interacting many- 
body systems is to some extent already determined by 
the rank of the interaction alone, and one does not need 
all the details of the interaction. We will show that odd- 
even staggering also fits into this picture and is not solely 
a consequence of an attractive pairing force. 

This paper is organized as follows. In Section || we 
introduce the Hamiltonian and discuss the odd-even ef- 
fects arising from the one-body part alone. Section |l| 
contains analytical results for the odd-even staggering 
due to a random two-body interaction (some technical 
details of this analytical analysis are included in an Ap- 
pendix). The crossover between the mean- field regime 
and the regime of strong interactions is numerically in- 
vestigated in Section The effects of breaking time- 
reversal symmetry are studied in Section |y| Finally, we 
give a summary. 

II. HAMILTONIAN AND STAGGERING 
INDICATOR 

As discussed in the introduction, we will consider en- 
sembles of Hamiltonians including only a single-particle 
energy and a random two-body interaction. We write 
this in the form 



M 

H = J2 £ i ( c It c 'T +4i c n) 

i=l 

+ C ^ Voaa'A^Aa* 

a, a' spin— pairs 

+ Cx ( 2 ) 

/3,/3' spin— 1 pairs 

The first term represents the mean-field contribution, 
where Si is the single-particle energy associated with 
orbital ?, and Cif, Cii are the one-particle annihilation 
operators for that orbital. As usual, we assume an or- 
dering £i < £j+i. The second and third terms repre- 
sent the interaction for pairs having spin S equal to 
and 1, respectively. The operators A a in the second 
term are spin-singlet two-particle annihilation operators 
A a = (ci[Cj^ — CiiCji )/y / 2(l + Sij) with a standing for 
the set of orbital pairs ij. The Ap in the third line are 
similarly defined for spin-triplet pairs. 

The randomness assumption tells us that there is no 
preferred basis within either the S = or S = 1 sector 
of two-body states. The couplings v saa > then should be 
taken from the Gaussian orthogonal random-matrix en- 
semble (GOE). We fix the variance of the v s to be unity 
for off-diagonal elements. The GOE then satisfies 

( v La>) = 1 +<W , (3) 

where (• • •) indicates an ensemble average and similarly 
for We are concentrating for now on the case of 

time-reversal symmetry, so the matrices vq and v\ are 
real and symmetric. The case of broken time-reversal 
symmetry in the presence of a magnetic field will be con- 
sidered in Section [v| 

The prefactors Co and C\ allow us to consider arbi- 
trary strengths of the spin-0 and spin-1 couplings relative 
to each other and relative to the single-particle level spac- 
ing. As we will see below, several qualitatively different 
regimes for ground-state staggering are possible within 
this simple random model, depending on the values Co 
and Ci as well as on particle density. 

Let us denote the ground-state energy of the A-body 
system as E{N). A useful staggering indicator is the 
empirical pairing gap 

A{N) = -[E{N +1)-2E{N)+E{N -1)}. (4) 

This three-point observable is essentially the "curvature" 
or second derivative of the binding energy with respect to 
particle number A. Positive (negative) A(A) indicates 
that the binding energy of the A-body system is larger 
(smaller) than the arithmetic mean of the binding ener- 
gies of its neighbors. We have an odd-even staggering 
whenever A(A) staggers with A. 

It is instructive to consider the trivial case where resid- 
ual interactions are negligible, i.e., Co = C\ = 0. Then 
the A— particle ground-state energy is given by E(N) = 
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2Et 7 i 2 e< for W even and E(N) = £(AT - 1) + £(jv+i)/2 
for N odd. Here TV may range between and 2M, where 
M is the number of available orbitals. One obtains for 
the empirical pairing gap 



A(N) = 



(£(jv/2)+i - s N / 2 )/2 > for N even, 
for N odd. 



(5) 



Thus, there is a trivial odd-even staggering due to the 
mean-field alone. In what follows we will mainly be in- 
terested in the effects of interactions, and in the effects 
of adding a magnetic field. For odd-number systems, a 
nonzero value of the empirical pairing gap must be due 
to interactions, and this allows one easily to discriminate 
mean-field effects from interactions. Such a discrimina- 
tion is more difficult for even-number systems and has 
recently been studied in mean-field plus pairing Hamil- 
tonians |i|,^,|3(J. We will see in Section fv| how mean- 



field effects can be distinguished from staggering caused 
by complex (or random) interactions. Note that an elec- 
tric charging energy -E c har g c = cN(N — 1) leads only to 
a N- independent constant shift A (AT) — > A(iV) + c and 
can therefore be neglected. 



III. EFFECTS FROM RANDOM TWO-BODY 
INTERACTIONS 

We now imagine the opposite situation from that of the 
previous section, i.e., we consider the regime Si = where 
mean-field effects are negligibly weak compared with the 
random two-body interaction. In this limit one might 
assume that all odd-even effects should disappear. Sur- 
prisingly, this turns out not to be the case. Instead, we 
find persistent odd-even staggering arising only from the 
random two-body interactions; stronger binding energies 
for even- TV systems are typically obtained in numerical 
simulations. 

To understand this result analytically, we first note 
that the spectral density of a system with two-body inter- 
actions approaches a Gaussian shape in the many-body 
limit TV —>■ oo |3l],|32|]. The ground-state energies for dif- 
ferent particle number or s pin sec tors are then largely 
determined by the widths VTrH 2 of the corresponding 
Gaussians, scaling as 



E « bVTrH 2 



(6) 



where it is assumed without loss of generality that TrH = 
0. The prefactor b depends of course on the details of the 
deviations of the spectral shape from an exact Gaussian 
form, since these deviations cut off the tails of the Gaus- 
sian. Following an analysis along the lines of Ref. 
where the spectral shape is expanded in terms of Hcr- 
mite polynomials, and then estimating the coefficients of 
these polynomials, one may conjecture that the prefac- 
tor b should scale as log N with the number of particles 
in the system. In any case, for our purposes it is suffi- 
cient that this prefactor varies smoothly with N without 



significant staggering, which is confirmed by numerical 
simulations. Eq. (H) is known to provide a good qualita- 
tive explanation for some observed behavior of low-lying 
spectra, even for moderate numbers of particles where 
the Gaussian approximation is far from valid. For exam- 
ple, a comparison of TrH 2 for different spin sectors helps 
to explain J = total spin dominance among the ground 
states of random interacting many-body systems |12],p7| ■ 



A. Dilute limit 

Applying this approach to the present problem, we 
need then to understand how TrH 2 depends on the num- 
ber of particles and other parameters of the system. For 
simplicity, we consider first the dilute limit A" <C M with 
a pure S = two-body coupling (Cj = 0). 

From previous work, it is known that for even A^ the 
ground state comes always from the sector of total spin 
J = 0. In the dilute limit, a typical basis state in this 
sector has the form 

N/2 

|*j=o) = 2-»/ 2 nCo^oJ^ - aJ, T oJ a )|0) , (7) 



where the A~ orbitals i z , j z are all distinct. One easily 
checks that the number of S = pairs in this state is 
(N 2 + 2N)/8, since the particles on orbitals i z and j z for 
a given z are in an S — combination by construction, 
while the remaining (N 2 — 2N)/2 pairs have a probabil- 
ity 1/4 of being in a singlet combination. Any of these 
5 = pairs, labeled by a' in Eq. (||), may be annihilated 
by the Co term in the Hamiltonian. Another S = pair, 
a, must then be created; there are M 2 /2 choices for a in 
the dilute limit. Thus, simply by counting the number of 
terms in the Cq part of the Hamiltonian in Eq. (J2J) that 
may act on a total spin J — basis state we find 



Tr( evcn n<s:m)H 2 — Cq 



M 2 (N 2 + 2N) 



(8) 



for Af even and Jf«M. 

Similarly, for odd A" the preferred many-body ground 
state has total spin J = 1/2. The typical basis state has 
the form 

N/2 

i*,= 1/2 ) = 2- w / 2 4 T n(«W.T - <A0\°) > 0) 



where we take J z = +1/2 without loss of generality, and 
the indices i z , j z , and k are all distinct. This state con- 
tains only (A^ 2 + 2Af — 3) singlet pairs, resulting in 

Tr 2 _ 2 M 2 (N 2 + 2N-3) 

tr(odd JV<Af)-n — O . (IV) 

The 0(1/N 2 ) difference in the widths explains the odd- 
even staggering in ground-state energies. Intuitively, the 
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result is easy to understand: the ground state of the odd- 
TV system is forced to have a slightly higher total spin, 
resulting in a slightly smaller fraction of spin-0 pairs and 
consequently a smaller effect of the Co term in the Hamil- 
tonian. This in the end is what leads to weaker binding 
for the odd-iV system. 

The above analysis also gives a quantitative prediction 
for the size of the staggering effect. Assuming in accor- 
dance with Eq. (||) that the ratio of ground-state energies 
is proportional to the ratio of the widths, we find 



\E 



even TV 



= \E , 



ddN\ 



1 



2N 2 



for large N in the dilute limit and therefore 
A(N) Co = (-l) N ^\E(N)\ 



(11) 



(12) 



to leading order. We may compare this with the size of 
the pairing gap for the mean-field dominated system. In 
the previous section, we saw that A(JV) = A/2 on aver- 
age for N even, where A is the mean level spacing of the 
single-particle spectrum. This can be normalized, how- 
ever, in units of the binding energy. This binding energy, 
i.e., half the many-body spectral width, is \E\ w MN A/2 
in the mean-field case. So the average pairing gap has the 
size 



A(iV) r 



i-ficld — 



MN 



\E(N)\ 



(13) 



for even N, which surprisingly is smaller than the pure 
interaction-induced pairing gap in the dilute limit N <C 
M. 

At finite particle density p = 2N/M, mean-ficld- 
induced and interaction-induced stagger are of compa- 
rable size, a characteristic difference being the vanishing 
of the pairing gap A(iV) for odd N in the mean-field 
case, Eq. ([5]), which is absent for the pure interacting 
theory. In addition, in the presence of fluctuations in the 
single-particle spectrum, mean-field induced A(iV) will 
itself fluctuate between successive even values of N, while 
interaction-induced stagger is predicted to be smooth. 
These analytic predictions will be verified numerically in 
Sec. [v] below. 



B. General results for finite density 

The above derivation, though strictly valid only in the 
dilute limit, in fact provides a correct intuitive expla- 
nation of the stagger at any density for a pure 5 = 
two-body interaction. Handling the 5=1 interaction 
requires more care, since the qualitative behavior will 
depend strongly on the density p. We therefore need the 
exact expressions for TriJ 2 in various particle number 
and spin sectors. These expressions may be straightfor- 
wardly, though perhaps rather tediously, obtained by ap- 
plying the original Hamiltonian, Eq. (^|), to various basis 
states and evaluating the norm. 



The full results are presented in the Appendix. There 
we find that for a pure singlet random interaction, the 
prediction of Eq. ([l^) for the size of the staggering, ob- 
tained above only in the dilute limit, is in fact confirmed 
as a lower bound for arbitrary densities in the many-body 
limit N — > oo: 



(-l) N A(N) Co > 



2N 2 



\E(N)\ 



(14) 



The situation is more complex for a pure triplet cou- 
pling (Co = 0), since here the ground state may be a 
state of either minimal or maximal spin. In this case 
we see using formulas given explicitly in the Appendix 
that a critical density p cr it exists below which there is no 
staggering, while above which interaction-induced stag- 
gering of order \E(N)\/N 2 appears, just as in the singlet 
case. As the singlet coupling is turned on, p cr it decreases, 
reaching at Co = C\. Thus, odd-even staggering with 
stronger binding for even- TV systems is predicted to be 
a very general consequence of random two-body inter- 
actions, present for pure-singlet and pure-triplet interac- 
tions as well as in the intermediate case. 



IV. CROSSOVER BETWEEN MEAN-FIELD 
REGIME AND STRONG TWO-BODY 
INTERACTIONS 

The analytical results of the previous sections were ob- 
tained for pure one-body or pure two-body interactions. 
In this section we will study the odd-even staggering for 
the full Hamiltonian (||) numerically. To this purpose we 
draw the random matrices vq and v\ in Eq. (|2|) from the 
GOE and compute the ground-state energies of Hamil- 
tonian of Eq. (||) for several particle numbers N. This 
procedure is repeated many times for each N to obtain 
ensemble-averaged values for the ground state energies 
E(N) and the empirical pairing gap defined in Eq. (Q). 
In what follows we set the number of single-particle or- 
bitals to M = 10, and obtain ensemble averages from 200 
runs. The largest matrices of the ensemble have dimen- 
sion 63504; their ground states are computed using the 
sparse matrix solver Arpack p3|. 

We have to assign values to the single-particle energies 
£j of the mean field and to the coupling constants Co and 
Cl of the two-body interactions. We assume a mean-field 
spectrum with level spacings £j+i — £j that are Wigner- 
distributed. This is consistent with the assumption that 
our quantum dot or metallic grain has irregular shape. 
To study the transition, we multiply the single-particle 
energies with a factor cost/9 and set the spin-0 coupling 
Ca{ip) = svmp. Here ip is in the range <p € [0,7r/2] and 
thus parameterizes the transition from the mean field to 
the regime of strong interactions. The spin-1 coupling 
Ci is set to zero. Figure ^ shows the empirical pairing 
gap (^) as a function of particle number N for parameter 
values <f> = 0,7r/12,7r/2. 
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FIG. 1. Empirical pairing gap as a function of particle 
number for parameter ip — (full line), ip = tv/12 (dashed 
line), and ip = n/2 (dotted line; graph scaled by a factor 1/2 
for display purposes) shows the transition from the mean-field 
regime to strong interactions in the spin-0 channel. 



We see from Figure [j] that the odd-even staggering 
persists throughout this transition. In the absence of the 
mean field (</> = 7r/2), the staggering decreases slowly 
with increasing N and then increases again very close to 
the maximal filling, when the number of holes becomes 
small and p approaches unity in Eq. (A4). Its envelope 
depends smoothly on N if only even or only odd val- 
ues of N are considered. These qualitative results are 
fully consistent with the analytical predictions obtained 
in Sec. Ill and in the Appendix. The absence of such 



a smooth envelope thus indicates that the staggering is 
instead dominated by mean-field effects, as in the ip = 
line in Figure |l|. Similar observations have been made for 
pairing-plus-quadrupole in Ref. 0|. Note that the ran- 
dom interactions drive the empirical pairing gap A(iV) 
to negative values for odd N; in this sense the staggering 
is more pronounced in the presence of interactions than 
in the mean-field regime. Note also that the magnitude 
of the staggering itself contains only little information 
since the transition from the noninteracting to the inter- 
acting Hamiltonian does not correspond to a transition 
in a physical system. 

We repeat these calculations in Fig. ^ for the case of 
vanishing spin-0 coupling, Co = 0, and set the spin-1 cou- 
pling to Ci(ip) = simp. Again, odd-even staggering per- 
sists throughout the transition. In the regime of strong 
interactions the magnitude of the empirical pairing gap 
increases with increasing N for even N. The situation is 
reversed for odd values of N. Leaving out very small sys- 
tems (N — 3), the envelopes for even and odd N are still 
smooth enough to discriminate mean-field effects from 
interaction-induced pairing. 
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FIG. 2. Empirical pairing gap as a function of particle 
number for parameter cp = (full line), cp — n/12 (dashed 
line), and p — tt/2 (dotted line; graph scaled by a factor 1/2 
for display purposes) shows the transition from the mean-field 
regime to strong interactions in the spin-1 channel. 

Finally, we consider the case of equally strong spin-0 
and spin-1 couplings and set Co(tp) = Ci(<p) = sin^>. 
Figure [| shows that this case is qualitatively similar 
to the case of pure spin-1 coupling, since triplet pairs 
outnumber singlet pairs by a 3:1 ratio in the large- TV 
limit. Again, the interaction-induced staggering exhibits 
a smooth envelope and can therefore clearly be distin- 
guished from mean-field effects. 
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FIG. 3. Empirical pairing gap as a function of particle 
number for parameter cp = (full line), ip — -n/12 (dashed 
line), and ip — tt/2 (dotted line; graph scaled by a factor 1/2 
for display purposes) shows the transition from the mean-field 
regime to strong interactions. 



V. MAGNETIC FIELD EFFECTS 

BCS-like pairing results from strong correlations be- 
tween fermions in time-reversed orbitals. Thus, these 
correlations can be destroyed by a sufficiently strong 



breaking of time-reversal symmetry. Examples of this 
well-known phenomenon are the breakdown of electronic 
superconductivity in the presence of sufficiently strong 
magnetic fields and the reduction of pairing correlations 
in rapidly rotating and deformed nuclei. In this section 
we want to study how breaking time-reversal symmetry 
affects the odd-even staggering in systems with a random 
two-body interaction. Having metallic grains in mind we 
thus consider the effect of a magnetic field. To be defi- 
nite, we take a uniform i?-field in the ^-direction. This 
leads to Zeeman splitting and adds the following one- 
body term to the Hamiltonian 



M 



H B = t iBj2 



(15) 



which also breaks rotational symmetry, i.e., only the pro- 
jection of the total spin J z remains conserved. Here, /i is 
an appropriate constant. A second effect consists of the 
modification of the random two-body interaction. Pro- 
vided the time-reversal symmetry breaking induces split- 
tings that are larger than the mean level spacing, the 
random matrices VQ aa / and vippi in the Hamiltonian (|^) 
have to be drawn from the Gaussian unitary ensemble 
(GUE). Accordingly, Eq. (|) for the S — matrix vo a a' 
and the corresponding formula for the S — 1 matrix Vipp* 
have to be replaced by 



(\vQaa' 



1. 



(16) 



This reduces the variance of the diagonal matrix elements 
by a factor of two when compared to the GOE. Consid- 
ering the random two-body interaction alone, this effect 
introduces only small corrections of order 1/A 2 to the 
results presented in the previous sections and in the Ap- 
pendix. 

Let us consider the trivial case where residual inter- 
actions can be neglected. The independent pairing gap 
then becomes 



A(A, B) 



2 ( £ N/2^ 



Z N/2) 



for N even, 
for N odd. 



(17) 



The odd-even staggering thus decreases with increasing 
magnetic field and disappears when the Zeeman split- 
ting 2\iB equals half the mean level spacing (£j+i — £j). 
Note that Eq. ( |l7j ) ceases to be applicable for stronger 
magnetic fields. In the limit of very large B-fields, the 
ground state becomes spin polarized (i.e., has maximal 
spin J = N/2) and any odd-even staggering disappears. 
Note also that a breaking of time-reversal symmetry leads 
to a positive pairing gap at odd N and can thereby easily 
be distinguished from the effects of interactions. 

We now include again the random two-body interac- 
tions and compute the empirical pairing gap as the mag- 
netic field is switched on. The number of single-particle 
orbitals is M = 6. At vanishing magnetic field we assume 
an equidistant mean-field spectrum with unit spacing. 



The two-body random interactions have fixed couplings 
Co = C\ = 1/10. We add the Zeeman Hamiltonian ( |15| ) 
to the system and increase the Zeeman splitting 2/iB 
from zero to its maximal value (ei+i — £4) /2. Simultane- 
ously, the variance of the imaginary part of the random 
matrix elements is increased from zero to one, being held 
proportional to the Zeeman splitting. Figure ^ shows 
that the odd-even staggering decreases with increasing 
Zeeman splitting. The remaining staggering is due to the 
interactions, which are relatively weak in this example; 
the transition from the GOE to the GUE in the random 
two-body matrix is very mild. For strong two-body in- 
teractions the odd-even staggering remains strong when 
time-reversal symmetry is broken. Thus, the breaking 
of time-reversal invariance has only mild effects on the 
ground-state structure in strongly interacting systems. 
This finding is consistent with a recent study of time- 
reversal symmetry breaking in the nuclear shell model 
with random two-body interactions p3j. 




FIG. 4. Empirical pairing gap as a function of par- 
ticle number for various strengths of the magnetic field: 
2^,B/(e i+ i - Ei) = 0,1/4,1/2 (full line, dotted line, dashed 
line) . 



VI. SUMMARY 



We have shown analytically and numerically that ran- 
dom two-body interactions cause an odd-even staggering 
in interacting few-fermion systems such as small metallic 
grains or quantum dots. Interactions tend to smooth out 
the odd- A and even-A dependence of the pairing gaps 
and can thereby be discriminated from the non-smooth 
mean-field staggering. As expected, the breaking of time 
reversal symmetry leads to a decrease of the odd-even 
staggering; this trend can however be countered by suffi- 
ciently strong two-body interactions. 
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APPENDIX: 

The derivation of interaction-induced staggering in 
Section III was obtained in the dilute limit TV -c M. 



For general values of N and M and couplings Co and C\ , 
a straightforward counting procedure results in the exact 
expressions 



Tr 



,7=0 



H 



r 2 



3C 2 



(>4 



A 2 (2M- A) 2 



^[C^(2M 2 + MN - N 2 ) 



N 
16 



-3C 2 (2M 2 - 7MN + 3A 2 )] 
[C$(6M + N- 2Nd(l - d)) 



-3C 2 (10M - 13A - 2Nd{l + d))} 



N 
16 



[8C 2 - 24C 2 (2 + d)] 



(Al) 



Tr J=1/2 ff 2 = 



r 2 



3C 2 



04 



N 2 (2M - N) 2 



^[C 2 (2M 2 



MN - N 2 ) 
3C 2 (2M 2 - 7MN + 3 A 2 )] 
I- 9MN ~ N 2 /2) 

15MA + 31A 2 /2 + 2N 2 d(l + d))} 



-[C 2 (-3M 2 
+3C 2 (M 2 - 

1 



l6 [C 2 (-9Af + llA-2Ad) 
3C 2 (9M - 31A - 8Nd)} 



73Cf] 



Tr 



J=N/2 



H 2 



(A2) 

[N 2 (M - N) 2 - N(M 2 - 5MN + 4A 2 ) 



+ N{7N - Ml) - AN] 



(A3) 



In each of the three above expressions for minimal and 
maximal spin states, terms are ordered by their relative 
importance in the many-body limit at finite density, i.e., 
in the limit N -> oo with p = N/2M = const. The 
leading term is 0(N 4 ) in the many-body limit, and this 
leading term is seen to be manifestly symmetric under 
particle-hole exchange N — > 2M—N for the minimal-spin 
states (of course the maximal spin states J = N/2 exist 



only for N < M). At subleading order, the symmetry 
is broken due to anticommutation relations between the 
creation and annihilation operators in Eq. (0). At both 
leading and first subleading order, TrH 2 is clearly iden- 
tical for the J = and J = 1/2 states, indicating that 
the staggering can occur only at 0(1/N 2 ), entirely con- 
sistent with our dilute analysis in Section III. It is also 
at this second subleading order that we first encounter 
the dimensionless quantity d, which we did not need to 
consider in the dilute approximation, d, taking values 
< d < 1, represents the fraction of particles in the 
basis state that live on doubly occupied orbitals. 

As discussed above, for a pure S = coupling {C\ 
vanishing) , ground states come always from the sector of 
minimal spin, and thus we are led to consider the quan- 
tity 



Tr S=0 H 2 



Tr S=l/2^ 2 



3 - 6p(l - p) - 8p 2 d(l - d) 



Tr 



S=l/2 



H 2 



A 2 (l- p) 2 



> 



2N 2 ' 



(A4) 



where terms of higher order in the 1/N expansion have 
been dropped, and the last inequality is easily checked for 
all possible values of filling fraction p and double occu- 
pancy fraction d. Thus, our original estimate, Eq. (|Lq), 
obtained using the dilute approximation, is confirmed as 
a lower bound to the amount of predicted pairing gap, 



(-1)" A(A0c o > 



2iV 2 



\E(N)\. 



(A5) 



The situation is more complex for a pure 5 = 1 
coupling (Co = 0), since here the ground state may 
be a state of either minimal or maximal spin, de pend- 
ing on t he density p. Comparing Eqs. (A1,A2) with 
Eq. (A3) at leading order in the many-body limit, we 
see easily that J = N/2 is preferred at very low den- 
sity, p < pcrit = (5 — 2v3)/13 w 0.118, but as density 
increases a transition should occur to ground states of 
minimal spin. The preference for maximal spin at low 
density is obvious, since high-spin states clearly maximize 
the fraction of particle pairs with aligned spins (S = 1 
instead of S = 0). On the other hand, the physical rea- 
son for the transition to minimal-spin ground states even 
with a pure S = 1 coupling for p > p cr %t is that at high 
enough density there are relatively few other high-spin 
states that a given high-spin state can couple to. 

The relevant result for our purpose here is that at low 
densities there is no predicted st agg er in the many-body 
limit, in accordance with Eq. ( |A3j ), but for p > p cr u 
minimal-spin states again become dominant . A calcula- 
tion completely analogous to the one in Eq. ([Al tells us 
that once again the pairing gap A is positive (negative) 
for even (odd) N and proportional to 1/JV 2 times the 
magnitude of the binding energy. Thus, 



A(A) Cl = 



(-l) N A(N) Cl > 



0.027 

lV2~ 



(p < p cnt ) 

\E(N)\ (p > p crit ) (A6) 
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The above analysis generalizes easily to the generic 
case where the two coefficients Co and C\ are both 
nonzero. For Co > C\, ground states are always expected 
to come from the minimal-spin sector, leading to positive 
pairing gap A proportional to 1/N 2 times the binding en- 
ergy. For C\ > Co, on the other hand, there will be a 
transition between no pairing gap at low density to pos- 
itive pairing gap at higher density, the critical density 
Per it approaching 0.118 for C\ » Co and approaching 
at C = Ci. 
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